Effects of relative orientation of the molecules on electron transport in molecular 

devices 



ZHOU Yan-Hong 1 , ZHENG Xiao-Hong 2 , XU Ying 1 , ZENG Zhao- Yang 1 Q and ZENG Zhi 2 
1 Department of Physics, Jangxi Normal University, Nanchang 330022, China 
2 Key Laboratory of Material Physics, Institute of Solid State Physics, 
Chinese Academy of Sciences, Hefei, Anhui 230031, China 
(Dated: May 29 2006) 

Effects of relative orientation of the molecules on electron transport in molecular devices are 
studied by non-equilibrium Green's function method based on density functional theory. In partic- 
ular, two molecular devices, with the planer Au7 and Ag3 clusters sandwiched between the Al(100) 
electrodes are studied. In each device, two typical configurations with the clusters parallel and 
vertical to the electrodes are considered. It is found that the relative orientation affects the trans- 
port properties of these two devices completely differently. In the Al(100)-Au7-Al(100) device, the 
conductance and the current of the parallel configuration are much larger than those in the vertical 
configuration, while in the Al(100)-Ag3-Al(100) device, an opposite conclusion is obtained. 



PACS numbers: 85.65.+h, 73.63. -b, 36.40.-c 



The electron transport properties of molecules and 
clusters have attracted tremendous interest in recent 
years since molecular devices made from special kinds of 
molecules or clusters may find important applications in 
future electronic circuits With the advantages of ex- 
perimental techniques, for example, scanning tunneling 
microscope (STM) and mechanically controllable break 
junction (MCBJ), the measurement of current through 
nanoscale systems is now allowed. An amount of inter- 
esting behaviors, such as highly nonlinear I-V character- 
istics, negative differential resistance (NDR) and electric 
switching behavior, are found in various systems such as 
organics^ carbon nanotubes^ DNAjS^ etc. A distinct 
difference between the molecular devices and traditional 
conductors is that the transport properties of the molec- 
ular devices depend not only on the characteristics of the 
functional molecules or clusters themselves, but also on 
many other factors, such as contact geometry and gate 
voltagetS, In particular, in experiment, it is very difficult 
to determine how the molecule or the cluster is put be- 
tween the electrodesiS In other words, it is hard to know 
the details about the contact geometry. In fact, in exper- 
iment, generally, measurements are taken for many times 
and the average is performed to obtain the conductance 
of a molecular conductor so that the details of the con- 
tact geometry can be taken into consideration as fully as 
possible. However, in theoretical simulation, it is easy to 
design all kinds of possible contact geometries, for exam- 
ple, the distance between the central cluster or molecule 
and the electrodes, the orientation of the molecule rela- 
tive to the electrodes, etc. 

In this Letter, we study the effects of the contact 
geometry on the transport properties by two model 
systems: Al(100)-Au 7 -Al(100) and Al(100)-Ag 3 -Al(100), 
where both Au 7 and Ag 3 are stable clusters with planer 
structures. In a general measurement, the distance be- 
tween the electrodes is fixed while the orientation of cen- 
tral molecule relative to the electrodes is undetermined. 
Therefore, in our simulation, we will only investigate the 
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FIG. 1: Model structure of a two-probe system with Au7 
cluster coupled to two Al(100) electrodes: (a) parallel config- 
uration and (b) vertical configuration. 
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FIG. 2: The conductance curves (a) and I-V curves (b) of 
the AU7 system: the solid line for the parallel configuration, 
and the dotted line for the vertical configuration. The Fermi 
energy is set to zero (the same for all the following figures) 



effects of the orientation of the cluster relative to the elec- 
trodes, with the distance between the electrodes fixed. 
Particularly, in each model system, two typical configu- 
rations, with the cluster plane parallel and vertical to the 
electrodes, are studied. 

Our calculations are performed using a recently de- 
veloped first-principles package TranSIESTA-C, which is 
based on the NEGF technique. The TranSIESTA-C, as 
is implemented in the well tested SIESTA method, is 
capable of fully self-consistently modeling the electrical 
properties of nanoscale devices that consist of an atomic 
scale system coupling with two semi-infinite electrodes 
as shown in Fig.l. Details of the method and relevant 
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FIG. 3: The conductance curve (a) and I-V curve (b) of the 
Ag3 system: the solid line for the parallel configuration, and 
the dotted line for the vertical configuration. 
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FIG. 4: The PDOS of the two configurations of the A117 sys- 
tems at bias 1.2V: (a) for the parallel configuration; (b) for 
vertical configuration. 



FIG. 5: (The transmission coefficient of the two configurations 
of Au7 as a function of energy at the applied bias voltage 
of 1.2V: (a) for the parallel configuration; (b) for vertical 
configuration. 
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FIG. 6: The PDOS of the two configurations of the Ag3 sys- 
tems at bias 1.2V: (a) for the parallel configuration; (b) for 
vertical configuration. 



references can be obtained elsewhere 10,13 . 

The model structure for our theoretical analysis is il- 
lustrated in Fig.l. The Au7 or Ag3 cluster is sandwiched 
between the two Al (100) electrodes with finite cross sec- 
tion. The electrode is chosen from the perfect Al crystal 
along the (100) direction, and the number of atoms in 
each atomic layer is arranged as 5,4,5,4 Four Al atomic 
layers (5, 4, 5, 4) are selected for the electrode cell. The 
noble metal cluster AU7 or Ag3 together with four sur- 
face atomic layers in the left electrode and three surface 
atomic layers in the right electrode is selected as the cen- 
tral scattering region, as indicated by the two vertical 
lines. These atomic layers in the central scattering re- 
gion are thick enough so that the perturbation effect from 
the cluster to the electrodes is limited to these layers, 
and therefore the remaining part of the electrodes can 
be treated as periodic systems. The distance between 
the electrodes is fixed at d = 7. 71a.M. for the AU7 case, 
and 5.85 a.u. for Ag3 case. The average nearest-neighbor 
distance of Au-Au bond is 2.72 a.u.; it's 2. 67a. u. for Ag- 
Ag bond 14 . In our calculation, the convergence criterion 
for the Hamiltonian, charge density, and band structure 
energy is 10 4 and the atomic cores are described by norm- 
conserving persudopotentials. 

The conductance curves of the AU7 cluster in two dif- 
ferent configurations are presented in Fig. 2(a). From the 
figure, it is found that the equilibrium conductance in 
the parallel case is much bigger than that in the vertical 
case. To investigate its non-equilibrium transport prop- 
erties, the I-V curves of the two configurations of AU7 
cluster are also shown in Fig. 2(b). Obviously, the current 
of the parallel configuration is larger than that of the ver- 
tical configuration over the given bias range. With the 
increase of the bias, the difference in the current of the 



two configurations is getting bigger and bigger, and the 
current of the parallel configuration is as nearly 4 times 
as that of the vertical one at 1.4V. While for Ag3 cluster, 
the equilibrium conductance (shown in Fig. 3(a) and the 
current (shown in Fig. 3(b) of the parallel configuration 
are smaller than those of the vertical one. 

The different effects of the relative orientation of the 
clusters on the transport properties in these two devices 
can be understood by studying the coupling between the 
electrode and the molecule via the projected density of 
states (PDOS) because it can give us information on 
how strongly the molecule couples with the electrodes 1 ^. 
Without loss of generality, in both devices, the PDOS 
of the two configurations under bias 1.2V as an example 
will be considered. 

The PDOS and the transmission coefficient as a func- 
tion of energy of the two configurations of AU7 cluster are 
shown in Figs. 4 and 5. It can be seen in Fig. 4 that near 
the Fermi energy, the PDOS of the parallel configuration 
is much larger than that of the vertical configuration. 
Corresponding to Fig. 4, the transmission coefficient of 
the parallel configuration also takes larger values near 
the Fermi energy than that of the vertical configuration. 
Then we conclude that in the parallel configuration the 
AU7 cluster has strong coupling with electrodes via Au- 
Al bond, which makes the electrons with an energy near 
the Fermi level transmit through the cluster easily, thus 
give rise to a very strong transmission and large conduc- 
tance. But no Au-Al bond is formed and two vacuum 
potential barriers exist between the electrodes and AU7 
in the vertical configuration. Therefore the coupling be- 
tween the cluster and the electrodes is rather weak, which 
results in much smaller conductance for the vertical con- 
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FIG. 7: The transmission coefficient of the two configurations 
of Ag3 as a function of energy at the applied bias voltage 
of 1.2V: (a) for the parallel configuration; (b) for vertical 
configuration. 

figuration. The current / is calculated by the equation 
1 = X l-v dET(E, V b ), where the transmission T(E, V) 

2 

is a function of the bias V and the incident electron en- 
ergy E. At the applied bias voltage of 1.2V, the transmis- 
sion coefficient of the vertical configuration takes much 
smaller values than that of the parallel configuration at 
the energy range [-0.6,0.6] eV. Therefore, the current 
of the parallel configuration has much larger value than 
that of the vertical configuration. 

The corresponding PDOS and transmission coefficient 
for the two configurations of the Ag 3 case at bias 1.2 V are 
shown in Fig. 6. In contrast, the PDOS of the parallel 
configuration is smaller than that of the vertical config- 
uration. The transmission coefficient of the parallel con- 
figuration also takes smaller values near the Fermi energy 
than that of the vertical configuration (see Fig. 7). Thus 
the current and the conductance of the parallel configura- 
tion are smaller than those of the vertical configuration 
in the Ag3 cluster. The reason is that, when rotating 
the Ag3 cluster from a parallel configuration to a vertical 
one, no vacuum barrier is formed between the electrodes 
and the cluster. Furthermore, the coupling for the verti- 
cal configuration gets stronger than that for the parallel 
configuration, since electrons can transmit through the 
three Ag atoms in a parallel way, not a series way now. 

The change of the coupling strength when the cluster 



is rotated can also be displayed by the charge transfer 
between the electrodes and the cluster4£*i£ For the AU7 
case, in the parallel configuration, due to the stronger 
coupling, there is a charge transfer of 0.92e from the AU7 
cluster to the electrodes, while in the vertical configu- 
ration, it is just 0.61e. In contrast, for the Ag3 case, 
0.87e charge is transferred from the cluster to the elec- 
trodes due to weaker coupling for the parallel configura- 
tion while 0.92e due to stronger coupling for the vertical 
configuration. 

In summary, a first-principles method is used to study 
the effects of the relative orientation of the molecules 
on the electron transport in molecular devices. Two 
molecular devices, with the planer AU7 and Ag3 clus- 
ter sandwiched between the Al(100) electrodes are in- 
vestigated. We find that the relative orientation affects 
the transport properties of these two devices completely 
differently. In the Al(100)-Au 7 -Al(100) device, the con- 
ductance and the current of the parallel configuration 
are much bigger than those in the vertical configuration, 
while in the Al(100)-Ag 3 -Al(100) device, it is just the 
opposite case. The different effects of the relative orien- 
tation of the clusters on the transport properties in the 
two molecular devices are attributed to the different cou- 
pling in the molecule/electrode contacts. In the vertical 
configuration of AU7, two vacuum barriers are formed be- 
tween the electrodes and the cluster, thus the coupling 
becomes much weaker when the Au 7 cluster is rotated 
from the parallel configuration to the vertical one. How- 
ever, for the Ag3 case, on one hand, no barrier is formed 
when the cluster is rotated from the parallel configura- 
tion to the vertical one; on the other hand, all the three 
atoms of the cluster play the role of a bridge in a parallel 
style and strengthen the coupling between the electrodes 
and the cluster. Consequently, the conductance and the 
current are also strengthened in the vertical configuration 
of the Ag3 case. 
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